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Abstract 

Using a qualitative analysis based on the Hamiltonian formalism 
and the orthonormal frame representation we investigate whether the 
chaotic behaviour which occurs close to the initial singularity is still 
present in the far future of Bianchi VIII models. We describe some 
features of the vacuum Bianchi VIII models at late times which might 
be relevant for studying the nature of the future asymptote of the gen- 
eral vacuum inhomogeneous solution to the Einstein field equations. 



1 Introduction 

Many studies have been made of the asymptotic states of spatially homoge- 
neous cosmological models. Studies of the early-time behaviour have focussed 
upon the possibility of chaotic behaviour in the most general cases 0, [|], 
[|], while investigations of the late-time asymptotics have focussed upon the 
relative generality of cosmological models which evolve to become at least as 
isotropic as the observed universe 0] . In this paper we are going to consider 
some aspects of the late-time evolution of one of the four most general classes 
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of spatially homogeneous universe, see [Q], 0, [|], and |]T8[ for related 
studies. The Bianchi type VIII universe has not attracted the same interest 
from cosmologists as some other homogeneous models because it does not 
contain an isotropic Friedmann universe as a special case. However, it can 
become arbitrarily close to isotropy. It is primarily of interest because of 
its status as one of the most general homogeneous universes and hence in 
some sense, yet to be made precise, may offer a leading order approximation 
to some part of the general cosmological solution of Einstein's equations at 
late times, far from the initial singularity. In fact, if space is topologically 
compact it becomes the most general of the ever-expanding Bianchi type 
universes 0- 

In what follows we will study certain aspects of the late-time behaviour of 
vacuum Bianchi type VIII models. In particular, we will try to determine 
whether the chaotic Mixmaster behaviour which takes place close to the 
initial singularity still takes place in the far future as t ^ +00. We will show 
how a qualitative analysis based on the Hamiltonian formalism can reveal 
some basic features of the future asymptote of vacuum Bianchi VIII models. 
These provide a basis for studying the nature of the late-time behaviour of 
the general solution to the vacuum Einstein field equations. 

It is also of interest to note that late-time cosmological behaviour in gen- 
eral relativity corresponds to early times in the so-called Pre Big Bang sce- 



nario suggested by symmetries of superstring theory. It has been argued [10 



that plane-wave space-times are the generic initial states for these models. 
So in this context the study of the asymptotic future behaviour of general 
vacuum space-times is also of interest. Also, recent work [|lT| has pointed 
out that inhomogeneous models of low-energy bosonic field equations for all 
superstring models and M-theory { D = 11 supergravity ) exhibit an oscilla- 
tory approach to the initial singularity, so it would be interesting as well to 
study how these oscillatory singularities transform in type VIII universes at 
late times using a similar qualitative analysis, see also |T2I . 



We will first introduce the field equations for the vacuum Bianchi type 
VIII model and then discuss the orthonormal frame approach and the 
Hamiltonian formalism to analyse the possible future asymptotic states. 
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2 The field equations for Bianchi type VIII 

In the metric approach, where the metric components are the basic vari- 
ables of the gravitational field, we can introduce group invariant and time- 
independent frame vectors Cq such that the line element is given by 

ds'^ = -dt^ + gabWW^ 

where W°- are time-independent 1-forms dual to the frame vectors Ca- In the 
case of vacuum Bianchi type VIII we can introduce three time-dependent 
expansion scale factors as 

gab = diag{a'^,b'^,c'^) 
and Einstein's field equations become 

[{a' + cr - b'] (2) 



abc 2a?}p'c^ 

[(a^ - H'f - c1 (3) 



(abc)- 1 2 ,2^2 4 



abc 2a^b'^c^ 

abc , 
- + - + - = 4 

abc 

Here, the dot denotes differentiation with respect to the synchronous time t. 
If we replace the synchronous time t by a variable rj defined by 

dt = {abc)drj 

then the field equations (l)-(4) become 

2{\na)" = {b'' + c^f -a^ (5) 

2{\^b)" = {a^ + c^f -b" (6) 
2{\nc)" = {a^-b''f-c^ (7) 
(Ina^)" + (In 62)" + (Inc')" = {\na')\\nb'')' + ihib'")' ilri^)' + i\nc^)' ilria^)' 

(8) 
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where ' denotes differentiation with respect to rj. 

These equations are invariant under the transformation t] —rj. Now, 
as is widely exploited in the literature, Q], 0, to a good approximation we 
have abc ~ t as t — 0, which implies 

rj oc ln(t) 

and the invariance of the field equations (5)- (8) under rj —rj corresponds 
to invariance under 

1 

^ 1 

in (l)-(4): an apparent duality transformation relating early and late-time 
behaviour. If this symmetry were exact it would imply that the presence 
of chaotic behaviour at early times {t 0) means that there must also be 
chaotic behaviour as t — +oo, as suggested in ref. |^. A crucial point in the 
above argument is the accuracy of the assumption that 

abc ~ t 

One can use the field equations (l)-(3) to test this assumption, by looking 
at the second derivative of the metric volume element abc. Taking the sum 
of the equations (l)-(3) we get 

2{ahc){ahc)- = (a' - b'^f + 2a^^ + 2c%'' + c\ (9) 

Since the right-hand side is the sum of four positive terms, if abc 7^ we have 

{abc)- > (10) 

This shows that the relation abc oc t cannot be exact and so the duality 
symmetry argument cannot be applied to relate late and early time behaviour 
for the vacuum Bianchi type VIII model. However, the inequality (10) will 
be useful in deducing the late-time behaviour as t ^ +00. There is another 
consequence of (10) worth pointing out : if abc is a power law 

abc ~ r 

then (10) implies n > 1. Since dt = {abc)drj, we would have 
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r] oc n 7^ 1 

1 — n 

Note that near n = 1 there is a structural instabihty with respect to the 
variation of t with respect to 77. When n > 1 we have t — > as 77 — ^ —00 and 
t — +00 when 77 ^ 0. However, this differs form the case where n = 1, where 
rj oc ln(t) (which is a good approximation close to the initial singularity) since 
here 77 — * +00 as t ^ +00. If n < 1, then we would have rj a.s t and 
rj —>■ +00 as t ^ +00. Tiny variations on either side of 7i = 1 would make 77 
a bad choice of variable for an early-time analysis of Mixmaster behaviour. 
This may be of importance for numerical studies of Mixmaster behaviour 



1^, which are notoriously difficult to carry out because of the slow time 
variation of the oscillations compared with the growth in volume. 

Let u"" be the 4-velocity of the Hubble flow, the Hubble parameter is 
defined by 

This parameter can be expressed as follows 



where l{t) corresponds to a geometric-mean expansion scale factor. The 
deceleration parameter q is defined by 

// 

q = — — 

The Hubble parameter then satisfies the following equation 

+ (11) 

Defining the a dimensionless shear parameter by 

V2 _ 
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with 0"^ = IcTafcCr'^^, where aab is the shear tensor, we have the following 
simple relation between the shear parameter and the deceleration parameter 
in vacuum: 

q = 2T?. (12) 

Equations (11) and (12) can now be used to express the Hubble parameter 
as 

For the isotropic (S = 0) Milne universe we have Hit) = t^^ so equation 
(13) expresses the departure from the Milne expansion rate for vacuum ho- 
mogeneous cosmologies. 

We specialise to the Bianchi VIII universe and calculate the Hubble 
scalar defined above for u"" = (1,0,0,0) using the Christoffel symbols for 
Bianchi type VIII models [13], so: 



Equation (13) then gives the following expression for abc : 

abc = exp(/ j^^f ) (15) 

with s{t) = J E^dt 

From the above equation we can get (abc)'', as 

For this to agree with the inequality (10) we require 

9-3(1 + 2^2) > 

which implies 

< 1 (17) 

As we will see this is consistent with bounds on S given by equation (28) of 
the next section. 
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3 The orthonormal frame formalism 



We use the orthonormal frame formahsm of Elhs and MacCallum |]T5| . The 
orthonormal frame Ca is chosen to be invariant under the action of the isom- 
etry group with cq = u, so that the commutation functions Y^^, defined by 



depend on the synchronous time t and are the basic variables of the gravita- 
tional field. With this choice of frame, the line-element takes the form 



where u"" are time dependent 1-forms dual to the vector fields Ca- 

The physical state of the models can be described using the components 
of the shear tensor and the curvature anisotropies Uab which occur in the 



expression of the commutation functions, see Wainwright and Ellis [|TB[ and 



[p!7[] , ||18|| . In the above frame we define 

Oab = diag {ail, (722,(^33) 

and 

Hab = diag{ni,n2,n3) 
and the linear combinations 

(7+ = -(0-22 + 0-33) 
^- = ^(^22-^33) 

n+ = -(^2 + ^3) 

-^<-"^' 

Now scale out the overall expansion of the universe by introducing the 
expansion-normalised variables 

S± = f (18) 
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iV„ = I (19) 

The cosmological model is then fully described by the dimensionless variables 
S_, Ni, N2, N3). The variables Sj- describe the anisotropy of the Hubble 
flow and the Na describe the spatial curvature of the group orbits. If we define 
a dimensionless curvature parameter by 

then we can deduce from the generalised Friedmann equation the relation 

1 = + X (20) 
Note that the curvature parameter K is related to the variables A^^ by 

K = ^[Ni^ + + ATg^ - 2{NiN2 + N2NS + N^Ni)] 

We can introduce a dimensionless time variable r defined by 

- = - (21) 
and related to the mean expansion scale factor, l{t), by 

I = loe^ (22) 

These variables yield the following simplified evolution equations for the vac- 
uum type VIII state variables : 

E±' = -(2-g)E±-^± (23) 
iVi' = (g - 4E+)iVi (24) 
N2 ={q + 2E+ + 2V3E_)N2 (25) 

N3 ={q + 2E+ - 2\/3S_)A^3 (26) 

where ' denotes derivatives with respect to r and the functions S± are given 
by: 
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^4 



6 



[(N2-N,f-N,{2N,-N2-N,)] 



S. 



2V3 



{N2-Ns)iN^-N2-Ns) 



This system, together with equation (20), is an autonomous system of 
first-order ordinary differential equations for which a monotone function can 
be found with respect to the time variable r. Different Bianchi types lead to 
different restrictions on the ranges of the variables Na- For Bianchi type VIII 
we require A'"! < 0, A^2 > 0, A^s > and the set defined by these constraints is 
an invariant set for the evolution equations. On this invariant Bianchi type 
VIII set there exists the following monotone increasing function : 

Z = {N.N^N.Y (27) 

As a consequence of eq. (20) the following bounds can placed on the state 



variables |jT6 



< < 1 
- 1 < S± < +1 
< A^i^ < 12 
< A^i(A^2 - N-s) < 12 
< -2A^i(A^2 + N3) < 12 
< (A^2 - Nsf < 12 



(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



Therefore, the invariant set 5* for Bianchi type VIII is unbounded in 
such way that A'2 + N3 can assume arbitrarily large values. Together with the 
existence of the monotone increasing function Z, given by equation (27), this 
means we do not expect Bianchi type VIII orbits to approach an equilibrium 
point as r ^ +00. 
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4 Hamiltonian Cosmology 



As discussed in section 2, in the metric approach, we can introduce a set of 
group-invariant and time-independent one-forms W"", so that with an arbi- 
trary time variable t' the metric has the form 

ds^ = -N{t')dt'^ + gab{t')W''W^ (34) 

In vacuum, the Bianchi type VIII metric is diagonal and we write the metric 
components in the form ||^ 

ga, = diag{e'^\e'P\e'^^) (35) 

with 

A = - 2/3+ 
/32 = + /3+ + 
/?3 = + - 

In the Hamiltonian approach the spatial Ricci curvature defines a potential 
function y(/3+,/3_). 

The variable is related to the overall length scale for the metric (34) and 
to the Hubble parameter by the relations 

/ = ke^° (36) 

H = (37) 

The evolution of the Bianchi type VIII universe is then described by the 
motion of a point in the mini-superspace (/3°, (3^ , (3~) evolving in the potential 
V{(3^,(3~). The variable (3^ can be chosen to be the time variable so that 
(3^ = 1 and = II~^. Equations (22) and (36) imply that we have (3^ = r. 

This Hamiltonian approach allows us to study the early and late-time 
asymptotic behaviour of Class A Bianchi-type universes qualitatively. For 
vacuum Bianchi type /, the potential function vanishes, V{I3^) = 0, so that 
the evolution is a straight line in /3+-space. For Bianchi type VIII the 
equipotential curves of V^(/3+) are approximately triangular in shape, 0], 
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with an infinite open channel which starts at the vertex on the (3'^ axis 
and becomes increasingly narrow as (3+ +00. As r — >■ +00 (evolution 
away from the initial singularity) the walls of the triangle contract. If the 
universe point moves initially in a straight line (corresponding to Kasner-like 
behaviour) it will bounce off a potential wall and then move along another 
straight line (a transition to another Kasner-like behaviour). The triangular 
shape of the potential allows the evolution to be chaotic, which corresponds 
to the Mixmaster behaviour and is familiar from studies of the evolution 
near the singularity (r —00). But, what happens as r — +00 ? The 
Hamiltonian picture implies that the triangular potential contracts, so this 
would enable the universe point to bounce off the walls in a chaotic sequence, 
ad infinitum. We shall argue however from the inequality (10), or the bound 
(28), and the fact that N2 + A^s can take arbitrarily large values, that this 
chaotic future evolution cannot occur as r ^ +00, contrary to the claim of 
Halpern, [Q. 

In order to show this, we introduce two further equations, relating the 



S± to the [3 variables. Since |T^ 



where S-t coincides with expression calculated in the orthonormal frame for- 
malism, using equation (12), we have the relation 

I = S2 = + s J = (^)2 + (^)2 (39) 

The Bianchi type VIII model exhibits Mixmaster behaviour towards the 
past and this chaotic evolution is reflected in the infinite sequence of Kasner 
era's and bounces off the triangular potential walls in the Hamiltonian picture 



1^. When bounces occur, the function 



dr dr 

exhibits a series of peaks (a bounce produces an abrupt change in or 
By equation (39) this implies that should exhibit a series of peaks as well. 
A graph of the function E is given in Wainwright and Ellis [0]. Now, if 
the triangular potential contracts as r — +00, the interval between bounces 
decreases as well. This means that time derivatives of /3± could diverge (the 
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peaks of (3±{t) would become increasingly narrow) and so the function (40) 
would be unbounded. This would imply, by equation (39), that would be 
unbounded, which contradicts the inequality (17) or (28). In this picture of 
evolution in an eternally contracting potential, one can imagine two solutions 
which avoid this problem. 

First, it could be that /?±(t) is of the form 

P^(r) = e-''^ oscihe''^) + c, 

with a, 6, c constants. In the above the osc function could have a periodic 
or chaotic oscillatory behaviour. Clearly this could give increasingly rapid 
chaotic oscillations of decreasing amplitude and using equation (38) one finds 
a bounded, eternally chaotic, oscillating S±. However, defining the function 
Po(^) by 

— = (2 - g)po 
dr 



one can show the following relations in the Hamiltonian approach |16| 



where Na = ^aNa, ( no sum ) with rii = C^23, "^2 = C'^si and 71,3 = C^i2. 
Using equations (41), (35) and (12), one finds, 

\NiN2N^\ ~ exp (6r - 6 j {I - T.^)dT) (42) 

We see that if S± oscillates eternally and satisfies (17) or (28), then (42) 
implies that |A^iA^2^3| increases as r — > +00 (see end of section 3). Now, 
equation (30) implies that A''! is bounded. Since (3± tend to constant values, 
the expression for (3i in (35) and equation (41) with a = 1 implies that 



is bounded, and this means that equations (35) and (41) with a = 2 and 
a = 3 will imply that N2 and A^3 are each bounded. Hence the Bianchi type 
VIII invariant set would be bounded and this clearly contradicts the fact 
that N2 + N3 diverges as r — > +00 because of equation (42). 
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Second, the universe point could slow down, but then we could end up 
reaching a point where /3± — constant and — 0. One can easily see that 
in this case equation (23) will imply that equation (20) will not be satisfied. 
Indeed, S_ = will imply that N2 = for Bianchi VIII , so equation 
(23) requires S_ = 0. Equations (24),..., (26), (41) and (35) imply that the 
Na become constants, but then (23) requires 5"+ = 0, otherwise S+ does not 
satisfy the bound (17) or (28). In order to have 5"+ = we must satisfy 
either A^i = or A'^i = A^2 = ^3- The first case implies that the curvature 
parameter K = and this is inconsistent with the Friedmann equation (20), 
and the second possibility implies that Ni = N2 = = Q because of the sign 
restrictions mentioned in the previous section for Bianchi type VIII. This 
also leads to i^' = and contradicts (20) as well. Note that this argument is 
not necessarily true in the non- vacuum situation (see appendix). 

So what are the alternative possibilities consistent with the Hamiltonian 
picture? The infinitely long channel along the (3^ axis offers a solution. The 
universe point could evolve eternally between the two walls of the infinite 
channel, for instance there could be an infinity of oscillations between the 
two walls of the channel. This evolution would correspond to a non-chaotic 
oscillation of the universe point about the axisymmetric solution (see below) 
characterised by 



= 

In this case, because the channel becomes increasingly narrow as /?+ — >■ 
+CXD, we expect /3_ to undergo increasingly rapid oscillations of decreasing 
amplitude as r — >■ +00. Similar behaviour of the universe point when inter- 
acting with a channel of the Bianchi IX potential can be found in Ryan and 
Shepley [0] for a so-called 'mixing bounce'. Equation (42) also implies that 



^2 + ^3 — ^ +00 as r ^ +00 and so the orbits in the state space escape to 
infinity. 

Non-chaotic oscillatory behaviour for Bianchi type VIII seems also to 
be observed in the behaviour of the Weyl curvature of other Bianchi type 
universes according to Coley |]19| and in Wainwright et al. In this 



interesting work the late-time behaviour of Bianchi type VIIq models with 
perfect fiuid source was investigated, (see also refs. |jl|, 0, |2^. If the 
universe contains matter with pressure p and density p, obeying an equation 
of state 
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p = (7 - l)p 



with 7 constant, then if 7 < | the universe isotropizes with respect to shear 
but not with respect to the Weyl curvature. So even if the shear decreases at 
late times, the anisotropic Wcyl curvature dominates the dynamics, produc- 
ing a form of self-similarity breaking at late times. If | < 7 < 2, then T? is 
oscillatory as r — > -|-oo, with /3_ a decreasing oscillatory function, while the 
Weyl scalar diverges. It was further argued that since an open subset of a 
more general class of models (i.e. Bianchi type VIII) will display features of 
Bianchi VIIq models, the non- self- similar behaviour of the Weyl curvature 
could occur at some stage during the evolution of Bianchi VIII models. 

Let us introduce the dimensionless electric and magnetic parts of the 
Weyl curvature 

Eab — ^abjp 
Hab — Habjp 

For all Bianchi class A models, in the standard frame E^b and Hab are diag- 
onal and we define: 

E+^^{E22 + E33) (43) 

E_ = ^{E22 - Ess) (44) 
and likewise for and By defining 

N+ = ^{N2 + Ns) (45) 

7V_ = i(iV2-iV3). (46) 

We have the following useful relations between the state variables introduced 
in the orthonormal frame approach and the electric and magnetic parts of 
the Weyl curvature : 
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E+ = ^+ + _ + S+ (47) 

E_ = J:_- 2S+S_ + S_ (48) 

H+ = -^NiJ:+- 3A^_E_ (49) 

H- = -3A^_S+ + ^(iVi - 4A^+)S_. (50) 
If N2 + N2, becomes arbitrarily large, then by equation (50), we can have 

lim H_ = +00 (51) 

and the Weyl scalar, given by 

W = El + El + Hi + H^_, 

can become arbitrarily large as r ^ +00, see Wainwright 0] for a discussion 
of the self-similarity breaking. Similar behaviour for the Weyl curvature was 
found in perfect fluid Bianchi type VIIq universes with equation of state 
parameter satisfying 1 < 7 < 2 



5 The axisymmetric (LRS) vacuum Bianchi 
type VIII model 

Our qualitative analysis shows that the late-time behaviour of vacuum Bianchi 
type VIII models will be non-chaotic. Evolution could occur which exhibits 
oscillations about the axisymmetric (LRS) vacuum Bianchi type VIII model 
with S_ = 0. The asymptotic late-time behaviour of this solution is known 



21[|. The Bianchi type VIII LRS model is described by the following re- 
strictions on the state variables, and these define an invariant set for the 
evolution equations (23)- (26) : 

A^i < 
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N2 = N3>0 
S_ = 

This state space is unbounded, since N2 can take arbitrarily large values. 
However, we can use the following bounds on N1N2 (see inequalities (28)- 
(33)): 

< -N1N2 < 3 



We introduce a new variable Z defined by Wainwright, |21 



Z^ = -N1N2 

It can then be shown that the system of evolution equations (23)-(26) admits 
the following equilibrium point : 

A^i = 
2 

This corresponds to a fiat plane-wave asymptote, which has W = 0. More 
precisely, it is the Bianchi type /// form of fiat space-time, with metric 

ds^ = -dt^ + e{dx^ + e^^dy^) + dz^. (52) 

The Bianchi type Vlh diagonal plane- wave metric with h = —1 reduces to 
this line element for a particular choice of parameters. 

Since the above equilibrium point solution, (52), is of LRS Bianchi type 
III with n'^a = 0, it is a particular case of the following general metric which 
includes Bianchi types /, VIq, VIII, V and Vlh with n°'a = [l^ 
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where and go constants and X{t), Y(t), Z{t), w{t) are functions of the 
synchronous time t, and 

w{t) = -2b J xy-^z-^dt 

where b is constant. The hne element (52) has, 

6 = 0,X = y = t, Z = 1, ao = go = 

Now since this hne element is a future attractor for the LRS vacuum Bianchi 
type VIII, and since (53) contains non-LRS solutions, it is interesting to 
perform a stability analysis of this line element for non-LRS perturbation 
modes. The perturbation technique for special cases of (53) has already 
been developed by Barrow and Sonoda p2|. In this approach one defines a 
new time variable as 

du dt 

as well as the new variables F = j^, G = and H = ^ so that the field 
equations for (53) in vacuum become 

F' = 2{ao' + qo') - ^ - F{G + H), (54) 

G' = 2{ao' + aogo) + - G{G + H), (55) 

H' = 2{ao^ - aogo) - H{G + H), (56) 
with the first integral constraints 

= 2F{G + H) + 2GH - (6ao' + 2qo') (57) 

and 

2Fao = qo{G-H) + G + H. (58) 

The unperturbed solution is the line element (52) with Fo = 1, Go = 1 
and Ho = 0. We require equations for the following small perturbations. 
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F 



Fo + f, 



G 



Go + 9, 
Ho + h. 



H 



Using the equations (54)- (58), one gets the following (the equation for F is 
eliminated using the constraint (58)) : 



The eigenvalues associated with this problem are Ai = —4 and A2 = 
— 1 and so the solution (52) is asymptotically stable under small non-LRS 
perturbations consistent with (53), thus supporting the prediction that the 
universe point will evolve about the axisymmetric solution with = in 
the narrow infinite channel along the /?_(_-axis. Taking into account the results 
of the previous section, this means we expect the future asymptotic state of 
the full non-LRS vacuum Bianchi type VIII to have 




The fact that S_ vanishes is a consequence of the following : since /3_ — > 0, 
the line element (34) becomes axisymmetric and so the field equations (5), (6) 
and (7) lead to the LRS solution for the scale factors, thus /?+ — > ^, however 
from equation (41) one finds that if S_ oscillates and does not decrease, 
then the curvature parameter K can be unbounded, which contradicts the 
generalised Friedmann equation (20). Also, although S_ decreases, it can be 
shown to undergo increasingly rapid oscillations : if is bounded, oscillatory 
and decreasing, then equations (23), (38) and (41) imply that /?_ ( and thus 
E_ ) undergoes increasingly rapid oscillations. If the universe point oscillates 
about the axisymmetric solution, then the metric shares a feature of some 
gravitational-wave space-times, namely that the solution oscillates about a 
given form of fiat space-time as t ^ +00. 



g = -4g - Ah 



(59) 



h' = -h 



(60) 
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6 Conclusion 



A qualitative analysis based on the Hamiltonian formalism and the orthonor- 
mal frame approach gives the following information about the future asymp- 
totic state of vacuum Bianchi type VIII models : 

a. The chaotic Mixmaster solution close to the initial singularity evolves 
at late times into a simpler non-chaotic solution, contrary to the claims of 
ref. 0, and corresponding to the fact that the universe point evolves along 
the infinite open channel of the Bianchi VIII potential and oscillates about 
the LRS Bianchi type VIII solution in such a way that the shear parameter 
exhibits increasingly rapid oscillations as r — *• +00 . 

6. The form of the line element (34) tends to the Bianchi /// form of flat 
space-time because — 0. It is interesting to compare this with asymptotic 
states of other vacuum Bianchi types : Bianchi types IV, Vlh, and Vllh tend 
to plane-wave space-times; Bianchi type VIIq tends to the special equilibrium 
point A'^i = 0, A''2 = A'^s = constant, with S4. = —1 and S_ = 0, which in the 
state space is a line emanating from a point corresponding to the Taub form 
of flat space; finally, Bianchi type V tends to the Milne form of flat space. 
So it appears that different forms of fiat and plane-wave space-times play 
an important role in the future state of these eternally expanding vacuum 
Bianchi models and one might expect that they occur in a leading-order 
approximation to part of the general inhomogeneous vacuum solution to the 
Einstein field equations at late times. 
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8 Appendix 

For a non-vacuum homogenous cosmology the Friedmann equation (20) be- 
comes 

1 = + K + ^] (61) 



19 



where fl is the density parameter defined by (SttG = 1) 



with p the energy density of a perfect fluid with equation of state p = (7— l)p. 
It appears that if < 7 < |, the LRS Bianchi VIII models admits the flat 
Priedmann model as future attractor i.e. T, — Ni — N2 — — and Q = 1, 
so that the argument developed in section 4 does not hold. This is simply the 
well known cosmic no hair attractor of power-law inflationary universes whose 
matter content violates the strong energy condition. In fact, for non- vacuum 
Bianchi type VIII, the analysis of section 4 can be adapted as follows when 
the strong energy condition holds. Note first that the deceleration parameter 
is now given by 

g = 2E2 + ^(37-2)Q 

and so we have 5>0if37 — 2>0. For all 7 satisfying this inequality, one 
can apply the analysis of section 4 and the universe point will again evolve 
in the infinite channel as r — > +00. 

We end this appendix by a remark about the stiff fluid case with 7 = 2. 
This case is of interest namely because the stiff fluid is equivalent to the 
presence of a massless scalar field. The evolution equation for Q is 

^ = [2g - (37 - 2)]Q (62) 

Now when 7 = 2, the field equations (1),..., (3) are identical for vacuum and 
non-vacuum cases. So, the inequality (17) is valid and this corresponds to 

g < 2 



and so equation (62) implies that Q 
tends to a vacuum solution. 
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